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Abstract 

We evaluate the differential cross section for the electroproduction of lepton pairs off a polarized 

nucleon target in the generalized Bjorken region to leading power accuracy in hard momentum. 

We discuss the importance of this process for phenomenology of generalized parton distributions. 
• ^ , 

^ A special attention is given to the sensitivity of physical observables, i.e., diverse asymmetries, to 
$—1 ' 

a I their dependence on scaling variables that allows to map directly the functional two-dimensional 
surface of generalized parton distributions. 
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1 Introduction 



A quantum mechanical system is determined by its wave function. Acquiring the latter from 
theoretical considerations, like solving the Schrodinger equation, or experimental measurements 
allows one to predict any physical observable of the system. The bulk of experimentally accessible 
quantities is sensitive only to the absolute value of the wave function and the phase of the latter 
is essentially unattainable. To circumvent the difficulty one has to measure correlations of wave 
functions, — or more generally the density matrix, — where the phase difference of wave functions 
can be probed. The interference of a test system with a reference source, possessing an a priori 
known characteristics, serves the purpose of reconstructing the missing phase of the wave function 
and thus one acquires complete information on the quantum mechanical system in question. 

The nucleon represents a relativistic multi-particle quantum system in a bound state whose 
dynamic is driven by strong interactions. It is the subject of intensive studies for several decades. 
Recently, it was realized that the most powerful theoretical tools in the analysis of the nucleon 
structure are one-quark and one-gluon correlations, dubbed the generalized parton distributions 
(GPDs) [H 121 El- They are analogous to a field-theoretical generalization of the phase-space 
Wigner quasi-probability function of non-relativistic quantum mechanics 0E1, — a specific Fourier 
transform of a density matrix alluded to above. A GPD depends on several kinematical variables 
and it has to be mapped as a function of all of them in order to predict such a fundamental quantity 
as the angular momentum of nucleon's constituents, which is given by their second moment 1^1 



There are essentially three experimentally feasible processes that have a clear theoretical un- 
derstanding and which can be used for direct measurements of GPDs: electroproduction of the 
photon eN e'N''y which is sensitive to the deeply virtual Compton scattering (DVCS) ampli- 
tude 13 El IHl El , photoproduction of a lepton pair — £iN' and electroproduction of a 
lepton pair eA^ e'N'ii III1II2]- However, among these only the latter provides a setup necessary 
for an independent measurement of a GPD as a function of both scaling variables ^ and t]. The 
former two reactions cannot entirely serve the purpose of testing the sum rule due to the reality 
of the final- or initial-state photons, respectively, which leads to the restriction ^ = ^rj. The last 
process is the most challenging from the experimental point of view due to small cross sections 
involved and requires high energy of the beam and high luminosity, on the one hand, and full 
exclusivity of the final state, on the other. 

The process eA^ — e'N'ii will be the subject of our present study, — a generalization to the 
full complexity of a short preliminary note [TT]. The outline of the paper is as follows. After 
presenting the general structure of the cross section, kinematical variables and addressing the 
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Figure 1: Subprocesses contributing to electroproduction of muon pairs. 

issue of reference frames in the next section, we turn to the discussion of the factorization of the 
process within perturbative QCD in terms of GPDs in section |3| In section HI we perform the 
computation of the cross section. We derive first a generating function for the squared of the virtual 
Compton scattering (VCS) amphtude and its interference with Bethe-Heitler (BH) amphtudes. It 
has a completely worked out leptonic part while the hadronic piece is expressed in terms of Dirac 
bilinears involving Compton and electromagnetic form factors. Next, we evaluate the products of 
bilinears for an unpolarized nucleon target, leaving the longitudinally and transversely polarized 
options for appendices. In section we give a detailed discussion of diverse asymmetries to be 
used for extraction of GPDs. Finally, we conclude. 

2 Kinematics 

The reaction we are dealing with consists of three interfering processes, depicted in Fig. ^ with 
implied crossed contributions. However, only one of them is sensitive to the one-particle corre- 
lations in the nucleon when at least one of the photon virtualities is large compared to a typical 
hadronic scale. It arises from the virtual Compton scattering amplitude, shown on the left hand 
side in Fig. ^ The other two amplitudes represent the Bethe-Heitler background. Presently, we 
discuss the production of a lepton pair of a different flavor compared to the one of the beam, i.e., 
the muon pair. The consideration of the electroproduction of electron pairs requires the addition 
of exchange contributions due to identity of the electrons in the final state, i.e., k' 

2.1 Phase space 

The generic form of the cross section of the exclusive electroproduction of lepton pairs off the 
nucleon, e{k)N{pi) e{k')N{p2)e{£^)i{£+), is 

da = ^—\r\^dLlPS^, (2) 
4pi ■ k 
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where T is a sum of the amphtude of the virtual Compton scattering and two Bethe-Heitler 
processes, T = Tycs + '^Hi + '7bh2 ? displayed in Fig. [T] The four-particle Lorentz invariant phase 
space 

rfLIPS4 = dM^gdUFSs d<Pu , (3) 

is factorized, by introducing the integration over the invariant mass of the lepton pair M^g, into 
Lorentz invariant phase-space factors for the production of a heavy timelike photon off a nucleon 

rfLIPSa = (27r)^5(^) (k + - k' - - g,) {pj - M^) {k'') ^^6^ {ql - M^,) , 

(4) 

and its subsequent decay into a lepton pair 

d^u = {f - m^) 5+ {{q, - l^f - r4) . (5) 

A simple calculation gives for them 

= d^Bdydi-A^)d4> ^Pdfk 

respectively, where the solid angle of the final state lepton in the ii center-of-mass frame is df]^ = 
sin OidOidifi. Here = (p2 — Pif is the momentum transfer in the t-channel and y = p ■ qi/p ■ k 
is the incoming lepton energy loss. The conventions for angles are obvious from Fig. El Here we 
introduced the Bjorken variable Xb and we will denote the virtualities of the space- and timelike 
photons as 

^B = ^-^, qi^Ml,, (7) 

^Pi ■ qi 

respectively. The nucleon mass effects are encoded into the variable e = 2x-bMj^/ Q and the final 
state lepton velocity in the ii center-of-mass frame reads 



Extracting the lepton charge from the amplitudes, one gets for the cross section, expressed in 
terms of experimentally measurable variables. 



da 



16(27r)3 



dxBdyd{-A^)d(j)dMl^df2e . (9) 



2.2 Reference frames 

Let us discuss appropriate reference frames to be used in analytical computations. To have a 
finite Fourier series of the cross section in terms of the azimuthal angles, we have just introduced 
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Figure 2: The kinematics of the lepton pair production in elastic electron-nucleon scattering. The 
coordinate system with the 2;-axis being counter-ahgned to the spacehke virtual photon is termed 
TRF-I, while the one with z' along the three-momentum of the timelike photon is named TRF-II. 



above the target rest frame in which the 2;-axis is directed in the counter-direction of motion of 
the spacelike virtual photon, called the target rest frame I (TRF-I), see Fig. |21 In this frame we 
obviously have 

pi = (M^, 0, 0, 0) , gi = 0, 0, -ql) , (10) 

where 



Q 

= — 

e 



Q 



For the remaining four-vectors we find: 

• The outgoing nucleon momentum p2 = (-^2,^2) has the components 



E2 = M, 



N 



2M, 



N 



\P2 



-A2(1-AV(4M2,)) 



and the scattering angle of the recoiled nucleon is 



COS 6' 



4:XbMn\p2\VTT^ 
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(12) 



(13) 



• The incoming electron four-momentum is 



k = {E, k"", 0, k^) = E{1, sin 0^, 0, cos 9e) , (14) 

with 

ye ^/IT^ ^ ' 



and as a consequence sin^g = eyl — y — + e^. 

• The four- vector of the timehke virtual photon reads 



q2 ^ (u;2, V UJ2) , (16) 

with 

^^ = f + ^' ^ ^ k1 = - , (17) 



£ 2Miv ' 

To evaluate scalar products which arise in contractions of the leptonic and hadronic tensors, 
it will be convenient to transform all four-vectors to a frame, termed TRF-II, where the z'-axis is 
directed along 52- The latter is achieved by rotating the TRF-I z-axis along the three- velocity v 
of the timelike photon by the scattering angle 9^ (which lies in the hadron scattering plane), 

cosO^ = , . 18 

^ 2Qu;2vVT+^ ^ ' 

In TRF-II, q2 — {0J2, 0, 0, U2v), and for the other momenta we get 

qi — {uji, qlsm.9^, 0, — g'^cos^^) , (19) 

and 

k — E {1, sinOe cos 9^ cos ip-y — cos 9^ sin 9^, — sin 9e sin c/?^, sin 9e sin 9^ cos </7-y -|- cos 9^ cos 9^) , 

(20) 

with (p^ — TT + (j). In TRF-II the vector pi is unchanged. In these formulas, the sinus of 9^ is given 

by 

4:Xb{1-xb) +e 



sin 9^ = ^ ^ J- (A2 - A2 . ) (A^ - A2 ) (21) 

in terms of the maximal and minimal momentum transfer in the ^-channel. 



:^2VTT^^ ((1 - xb)Q2 - xbM^Y - e-'Q'M^ , (22) 
with - (+) corresponding to A^i^ (A^^^). 
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A boost from the timelike photon rest frame to the TRF-II along the direction of motion of 
the photon with velocity v, see Eq. flT7|) . yields 

£_ = ( -002(1 + v(3 cos6i) , -M^i (3 sm6iCOS(pi, -M^^ /3sin 6'^ sin 9?^, -002 {v + (3 cos6i) \ , (23) 

\ ^ ^ / 

where Qn and <y9^ are the solid angles of i- in the ii center-of-mass frame alluded to above. The 
vector is simply deduced by a reflection, yj^ ^ + vr and Qn ix — Qi^ from (or equivalently 
by the substitution /3 — — /5 ). 

Using the explicit form of four-vectors in the TRF-II, we can readily compute the invariant 
products. Namely, the most nontrivial scalar products are 

^■A = -^^^YT7^{(^' + ^'^) (l-2^cosy., + ^) -A^(^l-XB(2-y) + ^)|,(24) 
A / , /i#2\ , as n i^sing^cost^^ 

't; Q^^-2xbM|+. 

Q^ + x^A^ 



+A -H " cos6'£ ^. (25) 



where 



1 



K ^ , (1 - ^ - ^'^V4) (4a:B(l - xb) + e^) (A^ - A^,J (A^ - A^J • (26) 

Eq. ()24j) reduces to the known expression from Ref. jHj for the real final-state photon Mfj = 0. 
Finally, since the expression for ■ A; is too complex to be presented here, we will give below its 
expanded form which is used in all practical calculations. To make the results look symmetric, we 
introduce the variable 



1 _ Pi ■ l+f/3cos6'^ l + cos6'£ 



(27) 



y Pi-q2 2 2 

which varies in the interval 1 < y < 00. 
2.3 Symmetric variables 

For our subsequent use, we introduce the symmetric combinations of momenta 

1 = ^{Qi + (I2) , P = Pi+P2, A=p2-pi = qi-q2, (28) 
and the invariants built from them 

q' = -Q\ e = ^, V=^, (29) 
p ■ q p ■ q 
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with the latter two being the generahzed Bjorken and skewness variables, respectively. These can 
be re-expressed in terms of experimental ones, discussed in the previous section, via the equations: 

Q' = 1[q'~m^, + ^^ , (30) 

for the inverse resolution scales, and 

for the scaling variables. Note that and ^ can take both positive (Q^ > M"^^) and negative (Q^ < 
M^^) values depending on the relative magnitude of spacelike and timelike photon virtualities, while 
1] < and ~ 1] > 0. To complete the set of formulas, the inverse transformations read 

and 

_ (e - V)Q' - eAV4 

(i_^)g2_^A2/2- ^^^^ 

The bulk of results given in subsequent sections will be presented in the symmetric variables. 
The following scalar products are needed in course of evaluations. Neglecting the nucleon mass 
corrections ~ M^/Q^, the above formulas ()22|) simplify considerably 

^min ~ ^-^V-iN 1 9 ; ^max ~ V 



1 9 * — max ^ ^ / -I ' 

and the transverse momentum transfer admits the form A^^ ^ (1 — 77^) (A^ — A^j^). The l/{p-q)- 
expansion of the scalar products (j21I), (ESI), and £_ ■ k, keeping the leading and sub-leading terms 
only, results into expressions 

k-A^-^j{l-2Kcos!f^), (34) 

^- ■ A ^ y| (1 + 2ircos¥?^) , (35) 

^- ■ ^ + v){i - + li^ - v){i - y) 



+aJ (1 — y){l — y){^'^ — rf) cos{lp^ — (pe) + 2riK cos(p^ + 2riK cos^pi > . (36) 



with 




\±^J^^EMEI). (37) 
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We note that in Eq. ()36|1 we replaced the original square root J {\ — y){l — y){\ — 'rf/i'^)Q^ in 



front of cos(v9^ — v?^) by crQ /^y (1 — y){l — y){^^ — rj'^), since Q is positive for the kinematics 
we are considering, a is +1. However, below in discussion of relations between the amplitudes we 
will make use of exchange symmetries that induce the interchange ^ —C, and —>■ in the 
underlying formulas. Under these substitutions, a changes the sign and takes the negative value 
— 1. Making use of the results, we have just derived, we are now in a position to proceed with 
the computation of amplitudes. Before doing this, let us discuss in the next section factorization 
properties and the structure of the virtual Compton scattering amplitude. 

3 Factorization of Compton amplitude 

The amplitude Tycs is the object of interest since it involves the off-forward Compton scattering 
amplitude on the nucleon T^^^. This tensor 



T,, = tj dSe''^-'{p2\T{j,{z/2)U-z/2)} \p^) (38) 

is expressed as a hadronic matrix element of the time-ordered product of the two quark electro- 
magnetic currents 

= J2 Qq'^qizht.'ipqiz) ■ (39) 

q 

Here we extracted the absolute value of the electron charge so that Qq is the quark's fractional 
charge. 

3.1 Structure of amplitudes in the generalized Bjorken limit 

Let us demonstrate the factorization of the VCS amplitude using the example of a cubic scalar 
model^. Although theorems of this sort already exist in the literature it is rather instructive 

to have a fresh inspection. We only need to analyze the singularity structure stemming from the 
denominators of particle propagators. Numerators present in QCD case will not be relevant for 
our present needs. Complications due to the gauge invariance in the realistic case will be treated 
in a straightforward manner. We will disregard throughout the contribution of crossed diagrams 
since they do not bring any new insights into the factorization property of the hadronic tensor. 
The tree- level Compton amplitude, see Fig. |Sl(a), reads 



^The absence of a stable vacuum state in this model is irrelevant for the demonstration of generic perturbative 
properties of scattering amplitudes. 
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Figure 3: Lowest order perturbative diagrams exhibiting possible short-distance regimes con- 
tributing to the asymptotics of Compton scattering amphtude. Thick hues correspond to the 
highly virtual propagators. 

where the small correction to the scaling contribution is e = (2M^ — A^/2)/(p ■ q). From here 
it is obvious that the expansion parameter is l/{p ■ q) = HQ"^, not as it is suggested by 
the definition ()38|). Actually, this scale approaches zero when the virtualities of incoming and 
outgoing photons are of the same magnitude. When p ■ q ^ oo and ^ is fixed, one recovers the 
scaling coefficient function l/{x — C, +^0) convoluted with the "perturbative" GPD which, to this 
order, is if(o)(x, ?7, A^) = 6{1 — x). 

The one-loop expression for the amplitude displayed in Fig. 01 (b) is 



(1) = ^9 



1 



J{El+iOa) 



(27r)4A;2(A; + pi)2(fc + p2)2(A; + pi + gi)2 (Att)^ Jo 
(Oa = 0-0;) where 

L \ a / a a ) a 

+ (a2 + as) 



(41) 



a 



(42) 



and a = Z)j=i c^j. Integrating in the vicinity of ai — > 0, which corresponds to the large virtuality 
of the corresponding line in the Feynman diagram, — the propagator between the photon vertices, 
— one gets, to leading order in l/{p ■ q), the contribution 

1 



rpSDl 
^(1) 



{p ■ q) "'-1 X — ^ + iO 
where we introduced the one- loop GPD, which absorbs mass singularities. 



(43) 



Hmix, V, A^) = ^ ri[da,5(x-l + {l-v)^ + il + v)^)'' 
(47r)^ Jo j-J-2 \ a a J 



ilEi+iOs 



(44) 



Here a = 02 + 03 + 04 and Ei = Ei[ai = 0] does not depend on large scales. Summing the tree and 
one-loop contributions, one gets the factorized expression of the form (|43|) with being replaced 
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(a) 



Figure 4: Short-distance regimes corresponding to loop diagrams in Fig. El 



hj H = /7(o) + + . . .. One can easily convince oneself that the aforementioned perturbative 
expansion of the GPD //(x,?], A^) arises from the light-cone operators matrix element: 



We refer to appendix ^ for the definition of the light-cone vectors. 

Other short- distance and infrared regimes lead to power-suppressed contributions compared 
to the leading one see Ref. 0. More particle attachments from the hadronic line to the 

hard propagator obviously lead to a stronger dumping of amplitudes, except for the effect of lon- 
gitudinally polarized gauge bosons in gauge field theories. The latter restore the gauge invariance 
of the naive GPD (j45|l . with the 0-scalars begin replaced by the quark fields ip. The analysis of 
the most general hand-bag diagram, i.e., with a single propagator between the photon vertices, 
does not bring new complications compared to the already discussed one-loop example and can 
be treated in an analogous manner. In QCD case, due to the fact that quarks are fermions, the 
numerators of their propagators cancel the power of p ■ g stemming from the denominator, so that 

Due to the timelike nature of the final-state photon virtuality, the quark- ant iquark pair can 
form an on-shell intermediate hadronic state before annihilating into the heavy photon. This 
happens when the hard momentum is re-routed around the photon vertex as demonstrated in the 
diagram 121 (c). This configuration can potentially generate leading asymptotic behavior when the 
photon virtuality is low. Let us demonstrate that it is actually suppressed. The two-loop diagram 
El (c) has the form 



.e"«-^MP2|0 0+, 0) (e-, 0+, 0) bi) . 



(45) 




(46) 



/ 



d^k2 1 



X 



[a/? + A(/3-/?i)]2 
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with a = Ylt=i Q^J! /5 = Sf=i A- The exponential reads 
E2 = a- 



ai + A (1 - 


P ) 


a + A (1 - 





a2 



PiPz 



+/53(i-^)m,V (47) 

In the short-distance regime, i.e., ai — >• 0, /3i ^ 0, where we define a = 02 + «3 + «4 and 
P = p2 + p3, we get, see Fig. H (b), 

= /■ /■ ^. ^'f-^-f^W-''"-^' . (48) 

{p-qyJ-i Jo {1 - u){x - ^ + lOy ' ^ ^ 
where was given above, since Ei = E2[ai = = 0], while 

nii){u, Ml^ = 77^ n d(^-J ^ - f ) {^^^3 (1 - I) M,2,| (49) 



(47 

is the first term in the perturbative expansion of the correlation function 

n{u,Mj^) = iq^.j d'ze"^-^ j ^e-^^"- (0|T {0(O_, 0+, O)0(O_, e+, 0), j(^)} |0) , (50) 

where the "electromagnetic" current is i{z) = |(/)^(z). Therefore, we observe that this short- 
distance regime is suppressed compared to the leading one, Eq. (jlHjl . 

In the QCD case, the suppression of contributions due to the hadronic component of the 
photon is much milder than in the scalar example and it is only {p ■ q)^^^"^ compared to the 
handbag diagram. The structure of the reduced amplitude is the same as in Eq. (j48|) . however 
with only one power of the hard-scattering coefficient l/(x — ^-|-iO) being involved and the vacuum 
correlator is of the form 

n^{u, q2) = {q2fiq2u - M'^i g^^^ nin{u, Mf^ 

= I J d'ze'^-^- J ^e-^^^"- (0|T {v^(0„, 0+, 0)7^V^(0_, e+, 0), j^(^)} |0) . (51) 

The correlation function can be saturated by the p-meson and reads 

r2^ _ "^P fpi-^) , ^ n \ r d^ 



niu,Mff) = ^— ^ ^-4^^ TT + ^uil-u) , (52) 
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Figure 5: The leading asymptotic region in the Compton scattering amphtude (a) and the 
factorization of longitudinal gluons into the path-ordered exponential (b). 

where (pp{u) is the p-meson distribution amplitude normalized according to Jq du(pp{u) = 1, while 
gp/{4:n) = 2.36 ± 0.18, nip = 770 MeV and Fp = 150 MeV are the p-meson decay constant, mass 
and width, respectively. The second term on the right-hand side comes from the perturbative 
contribution to the correlator, known to two-loop order ^3], a part of which is dual to the p-meson 
in the interval s G [0, Sq] and is absorbed there. The parameter sq ~ 0.8 GeV^ is the continuum 
threshold^. Due to divergence in the correlation function, one has to use a renormalized expression, 
stemming from the subtracted dispersion relation, Ur^u, MJ^) = n{u, M"^^ — n{u, 0). The gauge 
invariance of the light-ray operator involved in the correlation function is restored by means of 
Wilson lines as discussed below. 

From this result it is apparent that there is an extra imaginary part, besides the conventional 
s-channel discontinuity, in the VCS amplitude associated with the production of an on-shell in- 
termediate hadronic state by the quark-antiquark pair. Due to the current conservation, which 
implies f^2u{—i+) = 0, only one Lorentz structure contributes to the leptoproduction ampli- 
tude so that 

n,iu,q2)^u{i.hM-h) = -n{u,Ml)u{i.)^M-U) ■ (53) 

We should note, that due to the inequality, in general, of the skewness and the generalized 
Bjorken variable rj ^ =1=^, possible complications due to the singular structure of the hard coefficient 
function do not arise since the pole in the former is away from the "turning" point in GPDs when 
one of the struck partons has zero momentum fraction. Moreover, it was demonstrated that GPDs 
are continuous at this point, i.e., have no jumps, so that the region does not present a problem in 
the case r] = either ^ UHl US] • 

We finally conclude that the leading region in the amplitude is depicted in Fig. El (a), where 

in addition to considerations of the scalar theory, one is allowed to attach an infinite number of 

zero-twist longitudinally polarized gluons to the hard part. They do not induce power suppressed 

^The approximation of the continuum contribution by a step-function threshold in the spectral density causes 
a divergence in the real part at M'?, — sq, which is spurious. 
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contributions and factorize into the Wilson line along the trajectory of motion of the struck quark 
in the hard subprocess, see Refs. jTHl [T71 ITH] and JH] for the most recent and complete treatment 
of this issue. As a result the quark field in the naive definition of the GPD gets replaced by 

0+, 0) ^ Pexp (^tg dz_A4z_, 0_, 0)^ 0, 0) . (54) 

Due to the unitary cancellation of the eikonal lines beyond the photon absorption and emission 
points, the path-ordered exponential extends only between the quark fields, see Fig. El (b). 



3.2 Lorentz decomposition of Compton amplitude 

The gauge invariant decomposition of the hadronic tensor was found in Ref. by an explicit 
twist-three analysis at leading order of perturbation theory [see Refs. j^H 122 for independent 
developments along this line] . Presently, we concentrate on the twist-two sector at leading order in 
QCD coupling constant. As a consequence, we will not discuss gluonic contributions, in particular 
helicity-flip gluon effects which introduce new Lorentz structures into the hadronic tensor. The 
Compton scattering amplitude T^^/ admits the following Lorentz invariant decomposition in terms 
of Compton form factors (CFFs) 

+ 7r^^e\paPp Qa I 9 tie - ^^^^ I { 9ux - I rj, A^) . (55) 

2p-q \ p-q2j\ p-qi ) 

Here 5*^^23 = 1 Xo simplify notations we will set in what follows Vi = V and use the longitudinal- 
longitudinal helicity^ amplitude 

Vl = - Vi (56) 

analogous to the longitudinal structure function Fl of deeply inelastic scattering. 
The Compton form factors are decomposed into Dirac structures via 

V = h+H + e+£, Vl = h+T-CL + c+Sl , A=h+n + e+£. (57) 

with Dirac bilinears which read 



'2M, 

A 



A,, 

N 



hp = u{j>2)ipi5u{pi) , = ^^m(P2)75^(Pi) , (58) 



^Referring to incoming-outgoing virtual photons. 
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where u is the nucleoli bispinor normahzed as u{p)u{p) = 2M]^ and 75 = (^J)- 

As we demonstrated in the previous section, the inverse s-channel energy p ■ q sets the dis- 
tance between the quark fields in the Compton scattering amplitude. Therefore, for the QCD 
factorization to be applicable one has to impose the condition 

(p-g) >max{M^, lA^I}. (59) 

In this domain, the Compton form factors are factorized, as we have shown above, into calculable 
coefficient functions and GPDs via, cf. ()43p. 

in, £) (e, r^, A^) = f_ dx {Hi, E,) (x, r], A') , 

i=u,d,s,G ^ 

{n,£){^,V,^')= E I' dxC^^\x,0{H,,E,){x,v,^'), (60) 

i=u,d,s,G ^ 

where the sum runs over all parton species. Analogous factorized formula holds for Tii and Sl 
with the coefficient function C^~\x,^) being replaced by Cl{x,^). To zeroth-order in the QCD 
coupling constant only the quark coefficient functions i = q 

CW(a;,^) = 91 ^ 91 + 0{as), (61) 

give a nonvanishing contribution. The Compton form factor Vl does not enter in leading order. 
This is a consequence of a (generalized) Callan-Gross relation, V2 = ^Vi + 0{as). Thus, Cl{x,^) 
starts at next-to-leading order, i.e., 0{as), as we emphasized above. Perturbative corrections to 
C^^) and gluon contributions are available at one- loop order ^^ 1^ 1^ I77| 1^. 



4 Cross section 

This section is devoted to the calculation of the electropro duct ion cross section with the leading- 
twist Compton scattering amplitude analyzed above. In the next section we compute a generic 
form of the squared amplitude, where the hadronic part is left untouched and the leptonic sector 
is fully worked out. Then this intermediate result is used in section for the evaluation of the 
cross section for an unpolarized nucleon and is further used for deriving the results in appendix 
IHl when the polarization of the target is available as an option. 

4.1 Generating function 

The square of the total amplitude, T = Tycs + + ^H2) involves three essentially different 
contributions 

= |Tvcs|'+X+|Tbh,+Tbh,P, (62) 
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the square of the virtual Compton scattering amphtude, — bihnear in Compton form factors, — 
the square of the Bethe-Heitler processes, — independent on GPDs and expressed solely in terms 
of elastic form factors, — and, at last, the interference term 

1- = ^Cs'^^Hi + ^CS^Ha + ^CS^Hi + ^CS^Ha , (63) 

which is linear in Compton form factors. 

For the electron beam, the separate contributions to the total amplitude read 

_ 

Tvcs = u{£_)-ff,u{-£+) u{k')-f^u{k) T^^ , (64) 

_ _ ( I 1 \ 

^BHi = -YX2 ^i^-hM-h) uik') 7^- -7^ + 7i^ J,, , ^ 7m w(^) J>y > (65) 



^BHa = u{k')-ff,u{k) m(£_) + l y ^ ^ 7^ ) m(-^+) J,y , (66) 



corresponding to diagrams (a), (b), and (c) in Fig. [H respectively, including the crossed contri- 
butions in the latter two cases. The VCS tensor T^,^ was given previously in Eq. while the 
nucleon electromagnetic current is parametrized in terms of the Dirac and Pauli form factors 

J, = {P2\j,mpi) = F,iA')h^ + F,{A')e, , (67) 

using the bilinears from (j3H|) . The two amplitudes Tycs and ^BHa change the overall sign when one 
switches from the electron to the positron beam, and so do the interference terms involving them, 
while the Tbhi does not. Obviously, both BH amplitudes are related by the interchange of the 
momenta k' ^ £_ and k ^ —£+■ Moreover, we find that the VCS and the first BH amplitude are 
even under the interchange of the produced leptons in the pair, while the second BH amplitude 
is odd. This symmetry property in the timelike DVCS plays an analogous role as the charge 
asymmetry in the spacelike case |lUj . 

The evaluation of separate terms yields expressions which are represented as a Fourier sum of 
a few harmonics in the difference of the azimuthal angles ipi — = ipi — cj) — tt . Namely, we get 
for particular contributions to the total amplitude squared. 

• The VCS amplitude squared 

IT-vosI^ = Q^^_i.) t + cos (n(^, - 0)) . (68) 

has the following expansion coefficients 

= 1(2 - 2y + y'){2 - 2y + f) (VV^ + ^^t) + 4(1 _ _ ^)i!_!^v^Vi , (69) 

''^^ = -^^il-y)il-y)ie-v')i2 - y)(2 - y) (vvl + V^Vt) , (70) 
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2{1 - y){l - y) {VV^ -AA^) , (71) 
= iy(2 -y){2-2y + f) (VA^ + ^V^) , (72) 



2' 

= -|V(l-2/)(l-^)(^'-^'M2 - ^) {^lA^ + AVl) , (73) 
b^^^ = . (74) 

Note that contrary to the DVCS case, due to the virtuahty of both the incoming and outgoing 
photons, the Lorentz structure accompanying Vl does indeed contribute to the cross section and 
generates, e.g., the coefficient a^^^. 

The interference of the VCS amphtude with the BH ones will involve lepton propagators from 
the latter which will bring [conveniently rescaled] factors in the denominator 

{k' + Af = -27]p-qViik), {k- Af = -2r]p- qViik) , (75) 
{i+ + Af = 27]p-qV3{i^), {i^ + Af = 2r]p-qV4{i-). (76) 

The expressions are rather lengthy and are obtained by means of substitution of Eqs. ()24j) and 
(j26p into the left-hand side of the above definitions. We also introduce the following shorthand 
notations for the structures involving the nucleon matrix element of the quark electromagnetic 
current to make the formulas look as concise as possible, 

S^^v(k- -q,] .J^-^(k~ -q,] ■ A gi • Jt , (77) 



y J p 

S2 = V (i- - iq2] -J^-^ie-- iq2] ■ A g2 ■ J^ (78) 

V y J p-(i\ y J 

Til = e^,„paq^,k^ApJl , 7^2 = e^i.p^g/_i.Ap . (79) 

p ■ q p ■ q 

With these results at hand, we find, similarly to the previous analysis of |7vcsP; the interference 
contributions from the VCS and BH amplitudes. 
• The interference Tycs'^Hi- 



^VCsTJh, = ^2^2^2(,2_^2)g4A2p^(fe)p^(fe) E + A^^) COS (n(^, - 0)) , (80) 



where 



al = 4(1 -y){l- y) (^S^V - iTZ^A + 2^^^^SiV^ 
-(2 - 2y + y^) {2-2y + f) [^SiV + r/7^l^) 
-2^(l-y){2-y){2-m-r]) (7^2^ + ^^sV^ j , 
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a\ = 2a^il-y)il-y)ie-V') (82) 
X I (2 - (2 - y) (siV + n,A - ^S,Vl^ - A^^[r]S2V - ^TZ^A) | , 

al = -4(e - r/)(l - y)(l - y)(5iV + 7^l^) , (83) 
bl = -y{2 - y){2 ~ 2y + f){^S,A + vTZiV) - - r])il - y)y{2 - y)'JZ2{V + Vl) , (84) 

b\ = 2a^{l-y){l-y){e-v')y{^ - v) (^lA + U.V + ^U.Vl^ , (85) 
6^ = 0. (86) 
• The interference 7vcs^h2- 

with 

al = 4(1 - y)(l - y) (^vS2V - ^Tl^A - 2^^^52Vl^ 
+ (2-2y + y^) {2-2y + f) (^^aV + v^^^A 
-2|(1 - y)(2 - y)i2 -m + V) (niA + |5iVz.j , 

al = -2a^{l-y){l-y)ie-v') (89) 

X I (2 - y) (2 - y) (^52V - 7^2^ - ^S2Vl^ + V - ^n,A) | , 

«2 = 4(e + r/)(l - l/)(l - y)(S2V - 7^2^) , (90) 



,2 



6^ = y(2 - y)(2 - 2y + + ^7^2v) - 2(^ + 77)^(2 - y)(l - y)7^l (V + Vl) , (91) 

bl = -2a^il-y)il-y)ie-v')y{^ - y) (^52^ - Tl^V - ^7^2VL) , (92) 
bl = 0. (93) 

The unpolarized parts of the two interference terms must obey a symmetry relation, since the 
BH amphtudes and are related by the exchange k ^ — and k' ^ Obviously, 
qi ^ — g2 under it, while the Bose symmetry ensures the invariance of the Compton amplitude 
()64p with respect to this replacement. As we mentioned above, all the amplitudes have definite 
symmetry properties under the permutation of leptons in the pair £_ and, thus, we take the 

advantages of combining both transformations together, the above with k ^ — and k' ^ i^. 
From the definitions of the four-vectors (|2Uj) and (|23|) one can read off, after some algebra, the 
complete set of substitution rules 

Q^^Q\ (T^-a, A^A, ri^ri, y^y, ^Pi^^. (94) 
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Next, we remark that the product of the BH propagators (fTSj) and ([7^1) obeys the symmetry 
relation 



V^Va{Q\ a', V, y, y, ¥^e) = V{P2{Q\ A^ -e, y,y,(P = ^i) . (95) 

The prefactors in the interference terms ()8Up and ()87|) are even under the transformation ()94j) . 
Moreover, from the definitions (f77|) . (f75|) . and (f7^ we conclude also that Si ^ — 52;7^i ^ — 7^.2 
with (jHH). Taking all of our results together, we deduce that the Fourier coefficients satisfy the 
equalities 



2 _ _ 1 



where V and A are even and odd functions in ^, respectively. This property is a consequence of 
the definitions (j57p and (j6Uj) . where in the hard-scattering amplitude one has to replace ^ — iO 
by -i + 20. 

4.2 Angular dependence of the cross section 

After giving the generic expression for the total amplitude squared, we will elaborate the hadronic 
part and present the result as a Fourier expansion in terms of the azimuthal angles, (f) of the 
recoiled nucleon and (p^ of a lepton in the lepton pair. Before doing so, we write down the general 
angular decomposition of the squared amplitudes, which results from the Lorentz structure of 
the leptonic tensors contracted with the hadronic ones. The harmonics, appearing here, can be 
classified with respect to the underlying twist expansion of the hadronic tensor (|38|) and the Fourier 
coefficients are in one-to-one correspondence with helicity amplitudes, defined in the target rest 
frame. Extracting certain kinematical factors in order to match the normalization adopted for the 
leptoproduction cross section of a real photon in Ref. j9j, the square of the VCS amplitude and 
its interference with the BH amplitudes as well as the squared BH amplitudes admit the following 
expansion in these azimuthal angles 

l^^c^l' = ^4.2S^.2^ E cos(n0) + sr^S(^,) sin(n0)} , (97) 



y'y\v' - e)Q'A^ to { ^ ^1^2(0) 



2f (1 — rj)e^ ^ \ y f \ 



y 




(c^(0) cos{nipe) + s^(0) sm{nipe] 



PlP2(cpi) v^"'^'^) cos(nv9f) + s^2(0) siniriipe) 
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3 ~ 

± Yl ■nJ'lyT> i^ni'^i) cos(n0) + sl^{<^e) sm{n(f))] 



n=0 



Here the + (— ) sign stands for the electron (positron) beam and in the expansion we have used 
the relation between the azimuthal angles = vr + 0. In turn, the Fourier coefficients for the 
squared VCS term {i = VCS), the interference term with the ffist BH amplitude (fHHjl {i = 1), and 
the squared BH amplitude ()65|) {i = 11) are expanded up to the second order harmonics in the 
azimuthal angle ipe of the lepton pair 



2 



= H {c<^cos(mv9£) + cs^„sin(m(^£)} , 

m=0 
2 

^ni^e) = {^<mCos{m<^i) + ssl^sm{m^e)} . (100) 

m=0 

While for the interference of the second BH amplitude (jMjl with the VCS one [i = 2) and its 
square {i = 22), an analogous expansion is performed in terms of the azimuthal angle (p. The 
interference of both BH amplitudes {i = 12) is analogous to Eq. (jlOOp . however, it contains now 
a third order harmonics in the azimuthal angle (pe. The Fourier coefficients lineary depend on 
the polarization vector of the nucleon, see Eq. ()B.3|) . At the edge of the phase space the overall 
coefficient in the BH amplitude gets corrected according to 

(i^^jiziiiiL (101, 

We emphasize that 1 / (r/^ — C,"^) expressions in the squared VCS (j^Tfj) and interference (j^Hj) terms 
are corrected in analogous manner to ensure their correct behavior in the limits — > and 
M^, 0. 

Finally, we remark that all BH propagators, defined in Eqs. ()75|) and (fTBj) . are even functions 
in the azimuthal angle (p: 

= r,{27r - if) for ^ = {1,...,4}, (102) 

and, thus, even and odd harmonics can be clearly separated from each other. It is also worth to 
mention that V3Vi as a function of the lepton-pair solid angles ipi and 9^ satisfy the symmetry 
relation 

V^V^iOe, = VsV^iTT - 9e, y^e + 7c). (103) 

For later use, we mention that as a consequence of this symmetry the integration over dOg in a 
symmetric interval around the point 6e = tt/2 gives for any definite symmetric moment in the 
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following characteristic cos-Fourier expansions (for any number r) 

' [VsV.Y „=oi2... I oos{2nipe) J \ r(^^,) = r(7r - ^^,) j ' ^ ^ 



where ■& < ^. 



4.2.1 Squared virtual Compton amplitude 

At leading twist, it turns out that l^vcsP depends only on the harmonics cos (n{ip£ — 0)^ with 
n = 0, 1,2. Consequently, we find in this approximation, with the help of the addition theorem 
cos [n{Lp£ — (/))]= cos{rup£) cos(n0) + sin^rupi) sin(n0), the following relations between the Fourier 



coefficients: 



vcs ^ vcs 



^C^m^ — cs^^^ ~ , (105) 

S^nm — ^^nm — , 71 ^ m . 

The nonvanishing Fourier coefficients cc^^^ can be easily evaluated from the generic Eqs. ()fi9|) - ()74|l 
and products of Compton form factors. For the unpolarized target one defines 

where the functions C{J^,T*) depend on the set of Compton form factors. For this case we find 
at leading order in l/{p ■ q) 

A' 

4M2- 



Cvv,unp(-^, = (1 - V^)^* - rf{n£* + £n*) -\TT7^,+rf]££\ (107) 

cIS'unp(-^, n = (1 - ri^)nH* - ri\n£* + in*) - v'i^££* , m 



C^S'unp(-^,-^*) =0. (109) 

The spin-dependent results, including both longitudinally and transversely polarized target op- 
tions, are presented in the appendix IB. 11 

Note that in the (spacelike) DVCS limit, i.e., when one sets rj ~ — ^, we retrieve our previous 
result from Ref. [H] 

^DYCS DVCS ^VCS _|_ /^VCS 
'-'unp '-'VV.unp '-'^.A.unp ' 

and analogous relations for the polarized case, see Eq. ()B.10|) . In this way, we find for the 
unpolarized target the following nonvanishing Fourier coefficients in the twist-two sector 

= 2(2 -2y + y'){2 - 2y + f) {C^^%J:F, ^ + C^S'unpl^, n) 
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1 fi 

+^(1 - - m' - r/^)C^Snp(-^L, -^2) , (110) 



= jsj{i-y){i-y){e-v'){2 -y){2- y){c^Snp(-^' + ^^Snpl-^^, -^l} (m) 

ccSp = 8(1 - - ^) {C^gLp(-^,-^*) - CTA%n,i^,n} ■ (112) 

All other coefficients are expressed making use of Eq. (jlUSp . Note, however, that the tensor-gluon 
contribution induces further second order harmonics, which are not displayed here since they are 
suppressed by a power of a^. The Fourier coefficients for the polarized target are collected in the 
appendix IB.ll 

4.2.2 Interference of virtual Compton and Bethe-Heitler amplitudes 

In the leading-twist approximation the following general relations between the Fourier coefficients 
of the interference term are established 

ss™^^cc™^ and csfj ^ -scfj for {nm} = {12, 21, 32} , 

^^nm — '^'^nm — ^^nm — ^'^nm — ^ ^^^^ 7^ TTi ih 1 , (113) 

where n, m+l < 3 for INT= {1, 2}. For the unpolarized target, five nontrivial entries in the case of 
the unpolarized lepton beam appear, namely, cc™ynp, cc^^^p, cc^^^p, cc^^^p, and cc™^jjp5 which 
will be supplemented by three further Fourier coefficients in the polarized-beam case: cs^^np, 
sCio,unp5 ^iid sc™ynp, while sc^^^p ^ sc^^^p ^0. To find their explicit form we have evaluated 
the products of Dirac bilinears. Again, for the unpolarized nucleon target we get 

[■^.|v_4Q^(i^|5^-4c,„.,^(^), ,114) 
f^Av = 4,Q=ii^ ^'°Hcv..,pOT . (115) 



[JZ2} [yK si 



sm I 



[S2} yyi {yKsiYi^i] 

Ha = 4g^^i^ c,,..,(^) . (117) 

[7^2J yy^ [yK cos ipij 

We introduced here universal electric- and magnetic-like combinations of the form factors inter- 
twined with CFFs 

Cv.unp = Fi'H — ^ , C^.unp = —^{^1 + F2)'H . (118) 
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Then from Eqs. ()80|) - (jHm) and ()98|) the following nonzero Fourier coefficients for the first interfer- 
ence term are evaluated in a straightforward manner: 

CC^l,unp = -8^(1 - y)(2 - 2/)(2 - y)^3fJe|Cv,unp(^) + C^,unp(^) - ^Cv,unp(-^L)|ai9) 
CsJi,,np = -8Ai?y(l -y){2- y)^53m |Cv,unp(-^) + C^,unp(-^) + ^Cv,unp(-^L)| , (120) 

cclo,unp = Sir 3?e |(2 - 2y + y2)(2 - 2y + y^) (^Cv,nnA^) - C^,unp(^)) 

-8(1 -y){l- y)^^^Cv,unp(^L)| , (121) 

Cci2,,np = 16i^(l - y)(l - y)^3^e {Cv,unp(^) + C^,unp(^)} , (122) 



S'^lO.unp 



-S\Ky{2 -y){2~2y + f)^m |cv,unp(^) - ^C^,unp(^)| , (123) 

CC^l,,np = -8^(1 - y)(2 - y){2 - y)i-^SRe|Cv,unp(^) - C^,unp(-^) - ^Cv,unp(^L) |(124) 

SC^l,unp = 8A]rt/(l - Z/)(2 - y)^53m |Cv,unp(^) - C^,unp(^) + ^Cv,u„p(-^l)| , (125) 

CcL,unp = 16i^(l - y)(l - y)^3fJe {Cv,unp(-^) - C^,unp(^)} , (126) 

7] 

suplemented by ssjg.unp^ ssgi u^p) csgi ^np; ^'^d ssgg u^p) which arise from Eq. piHj) . The correspond- 
ing expressions for the second interference term follow from Eqs. ()87| ) - (jU^ and ()98|). For the 
unpolarized lepton beam they can be obtained from the symmetry under the exchange k — £_ 
and k' as discussed above in section 14.11 

/2 2 2 2 2\ _/l 1 1 1 l\ y^y (^07\ 
|CCq]^, CC]^Q, CC]^2) CC21, CC32 J^^^ — |CCq]^, CCj^Q, CC^25 CC21, CC32 J^^^ C^-? ' \^'^') 

Here one has to keep in mind that Cy and are even and odd functions in ,^ — iO, respectively, 
and K{^,y) = K{—C,,y = y). It turns out that the remaining coefficients for the polarized lepton 
beam satisfy the following symmetry relations 



{cSoi,sc^i} -J {cs^i,sc^i} 



2/ - 1 r 1 1 1 2 y - 1 1 

" SCin„„„ = -Wz SCi 



unp y 1 _ y L ^^' ^^J unp 4^-4 ' --iu,unp \j 1 - y 



10,unp ~ \ -I ^""10, unp 



_^ • (128) 

We note also that K{^,y) ^ -^(y — 1)/(1 — y)K{—^,y). The explicit expressions of the Fourier 
coefficients for the second interference term are given in appendix IB. 21 together with the general 
structure and results for the polarized nucleon. 

4.2.3 Squared Bethe-Heitler amplitude 

The expressions for the Fourier coefficients of the pure BH term (j99|) are extremely lengthy and, 
therefore, will not be displayed here in an analytical form. In the following, we merely limit 
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ourselves to leading terms in the asymptotic expansion as Q'^ —>■ oo. Namely 
ccJo,unp ^ -2^^ (l - I) (1 - 2/)|(2 - 2y + y'){2 - 2y + f ) (^1 + (129) 

c^olunp - 2l±^ (l - (l - (l - ^) (1 - I (2 - 2y + y^) (2 - 2y + (130) 



4M2 

ccllunp - 4i±^ fl - 1) ^^il-y)il-y) ie-v') (131) 



X (1 - ,)(2 - ,)(2 - ,~) { (l + 3i) (l - (f,^ - + Ji^(F, + 



cc 



13,unp 



4^ (i - f)' (i - %-) ^^ii-y)a-y)ie-v'), (132) 



A 



2 



11 ~ 


4^+^ 


20,unp ~ 


1—7] 


11 ~ 


<^ { 


22,unp ~ 


1-r/ V 



X (l-y)(2-y)(2-y)(^F,2-^F|j , 

- ^ + f ) - 3 - - - (^^' - ^^^') ' ^^33) 

C4^,unp - 8^ ' ^ " f ) " ?) " ^^'^^ " ^^^^^ 

and the rest are expressed via the already known coefficients 

ppii _ (^ ~ 0^ „„ii _ii ~ ppii -LoILilppii n-ic;^ 

'-'-24,unp ~ _|_ ^-j2 20,unp ! °°ll,unp ~ '^'-'ll.unp ~^ ^ _ ^^^13,unp ' \^'^'~') 

^^13,unp ~ "-^IS.unp 5 °^22,unp '-'-'22,unp 5 °°24,unp '-^24,unp • 

One has to realize that this expansion is only valid if one stays away from kinematical bound- 
aries, e.g., ?/ <^ 1 is required. The reason for this is that the leading terms vanish with (1 — y) 
and sub leading corrections become important as y ^ 1. Contrary to the DVCS case, it appears 
that no partial cancellation occurs between the numerator and the denominator in the BH ampli- 
tude squared, so that in general the Fourier decomposition goes as high as up to the forth-order 
harmonics. 

We note that the Fourier coefficients for the second BH-amplitude squared simply follow from 
the symmetry under k ^ — and k' ^ i^. Since we extracted one power of ^ in front of the 
squared BH amplitude (jHH), we obtain the substitution rule 



^^nm nm 



and ss„^ — ss^,„ 
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IZ'-^ , (136) 



while the remaining variables {ri,A'^,Q'^} are kept unchanged. For the interference term of the 
first and second BH amplitudes the Fourier coefficients are 



cc 



,12 

00,unp 



X 



1 — 77 77 \ 



7] 



2rf 



(137) 



2 

4M2 2 



l-rf 



cc, 



12 

02,unp 



8l±I f 1 _ iVl - CVi - y)(2 - y)(i - y)(2 -y)U- ^"'"^ 



1 — ?7 \^ 

1 + a 



7]^ 



A2 



X |(2-2|/ + y2)(2-2y + y2) 



f2 / A2. ^ 

S [ -1 "mm 



CC 



,12 

20,unp 



CC 



12 

22,unp 



l-T] \ TjJ \ Tj' 



\ I -| mm 



5" 



4M2^^ 



2 A^^2\ 

4M2 V ' 
(138) 

(139) 



+ i!±^(F, + F.)^ 



10- 



1-772 

2 



1 — 7?^ 



(i^i + F^y 



[l-y){2-y){l-y){2-y)[Fl 



4M2^V, 
(140) 



1 — rjf] 



(141) 



X (l-y)(2-y)(l-y)(2-y) 1 



A? 



2 ^^2^ 
4M2 ^ 



27^2 



1-7/2 



(i^l + i^2)' 



12 12 
S^ll.unp ~ CCii unp 



X {(2-2t/ + t/2)(2-2y + y2) + 8(l-y)(l-y)}(l- 



A2 



SS22,unp ~ CC22,unp 



(142) 



(143) 



Again, under the exchanges k ^ — £_ and k' both amplitudes are odd, i.e., Tbhi ^ —Tbh27 

and so their interference term is invariant. Thus, the Fourier coefficients must satisfy the relation 



Ccl2 = _ccl2 



y^y 



and ss 



12 



-ss. 



12 



y^y 



(144) 



which is the case for our result. 



5 Measurements of GPDs 



Having computed the cross section, we will discuss now the observables that can be used for the 
experimental exploration of GPDs. 
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5.1 Extraction of GPDs from electroproduction of electron pairs? 

The major complication in the experimental measurement of the process is a rather small mag- 
nitude of the cross section which is suppressed by two powers of the electromagnetic coupling 
constant ctem compared to a typical deeply inelastic event. The other obstacle is the contami- 
nation of the heavy-photon events by the background of meson production. The latter can be 
circumvented in a relatively straightforward manner by avoiding the regions of M'^^ close to meson- 
resonance thresholds. However, this certainly also restrict the phase space in the measurements 
of GPDs. For a general discussion of this issue we refer to Ref. j7j. We also note that a numerical 
estimate of this contamination can be done by means of Eqs. ()51|) - (j3^ . Indeed, the contribution 
of the p meson resonance to the beam spin asymmetry turns out to be small in a perturbative 
QCD estimate 

Before turning to the discussion of possible measurements of GPDs in appropriate observables 
and presenting more quantitative estimates, we have to emphasize that a clear study of GPDs can 
be done in experiments in which the tagged flavor of the lepton pairs differs from the one of the 
beam. In case they are the same, the results deduced so far have to be supplemented by contri- 
butions in which the final electrons are interchanged. Under this exchange, the outgoing electron 
momenta jump places, i.e., k' ^ i^. This obviously yields an essentially different dependence of 
the VCS amplitude on the external variables. To employ the process with identical leptons for 
the extraction of information on GPDs, one has to ensure that the momentum flow in the quark 
propagator in the handbag diagram remains large. Indeed, the scalar product 

2-y-ycos^r / A^\l 

p-q=-p-q + (145) 

that sets the scale in the exchanged VCS amplitude remains large, however, it is now timelike. 
Note that here and in the following we denote the kinematical variables that enter the exchanged 
amplitudes with a prime. The power-suppressed contributions depend on all kinematical variables, 
especially, on y, y and both azimuthal angels and (pi. Besides {2 — y)/y > cos9i, which is fulfilled 
by the usual kinematical restriction y < 1, no other kinematical cuts are required to ensure the 
applicability of perturbative QCD. Moreover, we find that r]' is given by r] in leading order 

while ^' receives a strong dependence on the leptonic variables: 

^, ^ ; 2cos.,-,(i + oos^,) _ 2yr^v^....., _ 

2 -?/(! + cos 6'^) 2-y[l + cos9i>) ■ qj 

We point out that if 9i approaches the edge of phase space, i.e., 9i {0,7r}, the absolute values 
of the scaling variables in Eq. (jl47j) become identical |^'| ~ |^|. The conclusions we draw from 
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our kinematical considerations are as follows. There are no crucial difficulties in the application 
of perturbative QCD as long as p ■ g is large, however, the analytical evaluation of observables and 
further studies are required to find an "optimal" method to deduce the (^, 77) shape of GPDs from 
measurements of the reaction e^p —>■ e^p e+e~. 

5.2 Mapping the surface of GPDs 

The most valuable information on GPDs can be accessed in observables that arise from the in- 
terference of the VCS and BH amplitudes, since they are proportional to linear combinations of 
the former. Such observables can be measured in single lepton- and hadron-spin asymmetries, in 
which the whole BH term squared drops out in the considered leading order in ttem, as well as in 
charge and angular asymmetries. 

In the former case one accesses the imaginary part of the VCS amplitude, where the contam- 
ination of the squared VCS amplitude is expected to be small. Thus, to leading order accuracy 
in the QCD coupling as, see Eqs. ()61|) and ()60p. one measures the GPDs in the exclusive region. 
Schematically, we write 



j^=n..£L 

where the imaginary part S^mjF for Ti. (7i) and S (S) is a (anti) symmetric function of ^ (however, 
always symmetric in 77) 



The kinematical factors }Cf{4>, can be read off from the Fourier coefficients, presented in section 
14.21 and appendix |B1 where the squared BH amplitude, the scaled propagators Vi and, at the edge 
of phase space, also the K and K factors should be exactly taken into account. In principle, a 
separation of different types of GPDs is partially feasible by a Fourier analysis or fully possible in 
case when single target-spin asymmetries are also available. 

For the other two asymmetries, mentioned above, one uses the fact that the ffist (second) BH 
amplitude is even (odd) with respect to the interchange of both the lepton charge of the beam, i.e., 
e~ ^ e^, and the lepton's momenta in the pair £~ ^ i.e., with respect to the simultaneous 
replacement 61 71 — 61 and —>■ ipi + tt. On the other hand, the VCS amplitude is odd 
under charge and even under angular exchanges. A combination of both asymmetries offers the 
possibility to access the real part of the VCS amplitude without the contamination from the BH 
background. We point out that a certain separation of different types of Compton form factors 
can be done again by means of a Fourier analysis of unpolarized measurements, while the complete 



da^ — da^ 
dcr^ + da^ 







(149) 
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separation requires double spin-flip experiments. Generically, these observables are analogous to 
the one in Eq. ()148|) . where S^nijF is replaced by the real part. To leading order in as, the real 
part of the Compton form factors is given by a "dispersion" relation: 

^eJ^iC, V, A^) = -PV dx -^'^mj^(x, r], A^) + . (150) 

71 J-1 q — X 

The value of this integral is also sensitive to the momentum fraction Ir^l < |a;|, which is not probed 
in single spin asymmetries. 

The unequal masses of the incoming and outgoing photons allow to probe GPDs away from 
the diagonal |^| = \ri\, the only kinematics which is accessible in DVCS. The skewness variable 1], 
given in Eq. (jHH), depends besides the photon energy ui = /2x-qM^ = yE (with respect to the 
target rest frame) mainly on the sum of both squared photon virtualities, i.e., + M'^^, while ^ 
essentially depends on their differences — M"^^. The boundaries of the (^, rj) region, probed in 
the process, are set by the following kinematical constraints: 

• The skewness parameter lies in the region 7]^^^ < r] < 0, where the lower bound comes from 
the kinematical condition |A^| > |A^j^|: 

- r/^in < v/-AV(4M2 - A2) . (151) 

• The upper and lower value of ^ is a consequence of the quasi-real limit of the space- or 
timelike photon 

-\V\<^<\V\- (152) 

• The minimal attained value of |^| = stems from the condition ~ M"^^. 

In Fig. El it can be seen that the area of the surface probed in the electroproduction of the 
lepton pair is quite extensive. The three contours displayed in this figure embrace the areas 
to be probed for different values of the electron beam energy E, the lepton energy loss y, and 
the t-channel momentum transfer A^: (i) solid contour corresponds to E = llGeV, y = 0.5, 
and A^ = -0.3 GeV^ (M^ < < lOGeV^), (ii) dashed contour corresponds to E = 25 GeV, 
y = 0.75, and A^ = -iGeV^ (-4A2 < < 20GeV^), (iii) dotted contour corresponds to 
E = 40 GeV, y = 0.9, and A^ = -3 GeV^ (-4A2 < < 35 QeV^). The corresponding values of 
the Bjorken variable are computed with the formula x-q = Q'^/{2M]\fEy). It is obvious that the 
higher the energy of the lepton beam, the higher A^ are allowed with observed applicability of the 
perturbative analysis of the Compton amplitude, and thus the higher values of |?7| are achieved. 

We have addressed above only the case > M^^ which probes ^ > component of GPDs. For 
the reversed inequality, one gets information on the region ^ < and probes patches of the two- 
dimensional surface analogous to the previous case. The positive mass of the final-state photon 
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Figure 6: The coverage of the GPD surface (the valence component of the u-quark distribution, 
as an example) with electron beams of different energies: (solid contour) E = llGeV, (dashed 
contour) E = 25 GeV, (dotted contour) E = 40 GeV, as described in the text. [For simplicity, we 
discarded in this plot the change of the GPD with for different kinematical settings.] 

allows to directly access only the exclusive, or distribution amplitude-like, component \r]\ > \^\ 
of the function. The inclusive, or parton distribution-like, component with |?7| < |^| requires 
spacelike virtuality for the outgoing photon which arises in two-photon exchange events in elastic 
electron-nucleon scattering. However, since the hadronic tensor (j38|) enters now via a loop integral, 
the single spin asymmetry measurements cannot be used for a direct extraction of GPDs. 



5.3 Compton form factors 

The magnitude of asymmetries depends on the relative strength of the BH amplitudes with respect 
to the VCS one. To get a rough idea of what happens we consider two limiting cases of the space- 
and timelike DVCS. In case of the production of a quasi- real final-state photon off the unpolarized 
proton target, the approximation of the amplitudes to leading power in for — A^^j^ ^ — A^ ^ 
4M^ < Q2 and ^ < 1 gives 

|r™|2 ~ y2 Q2 Fi(A2)2 ' ^ ' 

for M"^^ and large Q^. One expects that |7Y(^,^, A^)| behaves like and, thus, the ratio 



is essentially (1 — y)\A \/y Q . Obviously, the kinematical suppression by \A \/Q can be 
removed in the small-y region and, therefore, the DVCS cross section can be extracted in collider 
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Figure 7: The imaginary and real part of 7i(^, r/ = —0.2, = —0.2 GeV^) plotted versus ^ in the 
left and right panel, respectively, for different GPD models: = b^'^'^ = oo and B^"^^ = 5 GeV~^ 
(sohd, dashed) as well as for 6™^ = V'''^ = 1 and B''''^ = 9 GeV"^ (dash-dotted, dotted) with 
(dashed, dotted) and without (solid, dotted) D-term. 

experiments. On the other hand, for moderate values of y, one finds a rather sizeable (beam-spin) 
asymmetries in fixed target experiments. In contrast, in case of the timelike DVCS, replacing the 
variables y ^ y and above equation, we find the ratio 



^DVCS 


|2 sin^^^, 






^BH 


r\j - 

2 4 




Fi(A2)2 



(154) 



for — and large M^^. It is always suppressed by l/iWj^ and so the BH cross section overwhelms 
the VCS ones JU] ■ Thus, one anticipates that the second BH amplitude dominates the cross section 
except when one approaches the limit M^i —>■ 0. On this basis on would naively argue that the 
asymmetries will be suppressed as well. To get rid partly of this problem, instead of integrating 
over the full range of the lepton-pair scattering angles, one rather has to sum over a restricted 
domain and exclude the endpoint regions 9e = {0, vr}. 

To give quantitative estimates and provide some insights into the procedure of extracting GPDs 
from such experiments, we use several GPD models with quite distinct behavior in the central 
region. For the GPDs Hg, where q labels the three light (anti-)quarks, we take the factorized 
version of the {C,,ri) and A^ dependence, i.e., Hg{C^,ri, A'^) = Fg{A'^)Hg{C,,ri, A"^ = 0), where the 
reduced GPD at zero momentum A^ = is modeled via a factorized double distribution (DD) 
ansatz ^H]- For specific details we refer the reader to Ref. [H]. The essential freedom left within 
this specification concerns the parameter b'^ of the profile function, which controls the strength of 
the ?7-dependence of the reduced GPDs Hg{^, r], A^ = 0), the slope of the partonic sea form factor. 




(155) 
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and the D-term^, given as an expansion in terms of Gegenbauer polynomials [HI] 

D{^,V,A^) = - l^'^^g"^^^^ . (1 - x^) {-ACrix) - 1.2Cf (.) - OACrix)) \ , . (156) 

- AV0.77GeV2) ^ ^ '==^ 



The latter is entirely concentrated in the central region |^| < \ri\. The imaginary and real parts 
of the Compton form factors for these models are shown in Fig. In case the skewness effect is 
eliminated (as b oo), we model the reduced GPDs by the usual forward parton densities taken 
in MRS A' parameterization at the input scale = 4 GeV^. Consequently, for this so-called 
FPD-model the imaginary and real parts strongly increase as ^ gets smaller. This is displayed 
by the solid and dashed lines in Fig. [7| On the contrary, for 6 = 1, — for brevity we call this 
choice as the DD-model, — (dash-dotted and dotted lines), the contribution in the central region 
is suppressed and goes to zero for the imaginary part when ^ ^ 0. At the same time, the real 
part approaches a constant, with the value determined by the inverse moment of GPDs 

siten{^ = 0,V,A') = -f- Y: QI{h,{x,v,A')-H,{-x,v,A')} . (157) 

^ q=u,d,s 

The amplitudes supplemented by the D-term contributions are presented by the dashed and 
dotted lines in Fig. [7| It is clearly demonstrated that their effect is especially prominent (within 
the parameter range chosen for the estimates) in the real part of the Compton form factors where 
it changes their sign in the vicinity of ^ —>■ ±77. 



5.4 Single-spin asymmetries 

Now, having discussed the properties of the Compton form factors, we will have a closer look on 
single-spin asymmetries, in particularly, the beam-spin asymmetry for the proton 

da^ - da^ oc (±r*^^ + Tb^h,) SmTvcs + ■ ■ ■ • (158) 

Potentially, the interference term could be contaminated by the imaginary part Q^mTvcs^cS' 
arising from the interference of twist-two and -three Compton form factors. As for the DVCS 
process Q we expect, assuming the smallness of three-particle correlations, that this contri- 
bution can be safely neglected. As mentioned above, in leading order of perturbation the- 
ory, the single spin asymmetries are directly proportional to the linear combination of GPDs 
T.q=uAs Ql T Fq{i, -r], A^)}, where (+) - applies for (axial-) vector- type GPDs. For 

instance, eight leading-twist observables are measurable in the beam-spin asymmetry, which are 

^Note that this D-teim, taken from the quark sohton model |29[ I3U| at a low scale fi ~ 0.6 GeV, will mix mider 
scale evolution with a gluonic D-tcrm. For the present considerations it is suffice to neglect the scale dependence 
completely. 
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coming in pairs from the interference of the VCS with the first and second BH amphtudes^: cSq^^, 
sc}o, SC21, csjg as well as cSqi, scfg, SC21, csfg, see Eq. fll28|) . However, they depend only on two 
different linear combinations (111811 of GPDs: 



F,{n + Ul) - ^F2{S + Sl) , -v{Fi + F,)n + I (^F.Hl - ^F2Si^ ■ (159) 

Consequently, there there exist six constraints among the whole set of coefficients, which can be 
expressed as: 



(2-y)i2-2y + f)^il-y){C-v) j , , 



SC21 scio csqi y^(l — — T]) 



(160) 



supplemented by the relation (jllHj) . i.e., CS21 — — SC21 and cs^i — — sc^i. Another consequence of 
Eq. ()159p is that the beam-spin asymmetry gives us no handle on the Callan-Gross relation, i.e., 
the longitudinal Compton form factors can not be separated from the leading ones. The Fourier 
coefficients are projected out by taking the following moments, when integrated over the scattering 
and azimuthal angles: 

sin (f) — > sc}q , cos 9i sin ipi cSq^^ , cos 9i cos ipe sin(20) scl^ , cos 6'^ sin (pi cos(20) — csl^ , 
sin (f)e — > sc^Q , cos 9i sin (p — > cSq^^ , cos sm{2(p>i) cos — > scgi , cos 9i cos(2y9^) sin csl^ , 

where the weight in the first (second) row is even (odd) under the refiection ^ n — 6i and 
ipi 7c + ifi. In the same line of thinking, one can study the Fourier coefficients for the single 
target-spin asymmetries. Because of the substitution rules ()B.34|) and ()B.35|) . the number of 
the Compton form factors will be the same as for the case of charge and angular asymmetries, 
discussed below. Of course, single target-spin asymmetries are given by the imaginary part of new 
linear combinations in GPDs with a characteristic angular dependence. 

As we conclude from Eq. p60|) . the same information on GPDs is obtained by taking the 
appropriate moments in or (pg. However, the size of the complementary beam-spin asymme- 
tries can vary. Moreover, if one takes the asymmetry from the interference with the second BH 
amplitude, the weight must be odd and, thus, we have no contamination from the squared VCS 
amplitude. To suppress the squared BH amplitude, we integrate over the region 7r/4 < 6g < 37r/4 
and form alternatively the sin0 or simpi moments. This picks up the coefficients sc\q and sc^g^ 



^For brevity we neglect throughout this section in the Fourier coefficients the subscript "unp" , which refers to 
the unpolarized target. 
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respectively, cf. Eqs. (j^Hj) . (jlOOj) . and ()113|) . Thus, the beam-spin asymmetry 

f At^* 1 1 /•3V4 _ i-2n ^ /-a^ , , f 2sin(/) 1 da^ - da^ 



JO Jo 



l^^sinvjfj J\f Jit/a Jo Jo [2sin(y9£j dQidcf) 

oc S3m - ^i^2^: + + F^)n^ , (161) 

with the normahzation factor being 

r3^/4 /.27r da^ + dcj^ 



M= deJ d^e / 

Jir/i Jo Jo 



In/ 4. Jo Jo dQidcj) 

is analogous to the one defined in the case of space- and timelike DVCS. For the proton target it 
is mainly sensitive to the contribution Fi^mH and, therefore, we might neglect in our estimate 
the other two Compton form factors. 

In Fig. IHl we show the beam-spin asymmetries p61|) for an 11 GeV electron beam and r/ = 
-0.2, A2 = -0.2 GeV^ (A^i„ ^ -0.15 GeV^), and y = 0.5. We fix the value of + ^ 
3.4 GeV^, i.e., Q^/$, ~ 8.5 GeV^, and scan the —0.2 < ^ < 0.2 region by varying the virtuality 
Q2 in the range < < 3.4 GeV^. In the spacelike DVCS limit (the left panel with ^ -t], 
xb ~ 0.33), we uncover a typical beam-spin asymmetry of order 20% measured in fixed target 
experiments. This asymmetry is getting smaller in the DD-model (dashed and dash-dotted line) 
with increasing ^, reflecting the fact that the imaginary part of the Compton form factor goes 
to zero (see Fig. [7j). We stress that this feature truly arises from the GPD model, since C,/Q^ 
is fixed and so no essential kinematical suppression of the interference term arises as ^ — > 0. 
The negative D-term slightly changes the normalization away from the edges of phase spaces. 
In case of the FPD-model, the interference term increases at ^ with ^^^^^ (for MRS A' 
parametrization A = 0.17), as a consequence of the small-x behavior of the sea-quark densities. 
On the other hand, the VCS cross section in the denominator increases like ■C"^"^^? overwhelming 
the BH contributions, and forces the asymmetry to vanish at ^ = 0. This explains why with the 
FPD-model the asymmetry remains sizeable over a large interval of ^. The asymmetry, formed 
with sin0 is considerably smaller in the timelike region, compared with the spacelike one, since 
the first BH amplitude, responsible for this asymmetry, is getting smaller, while the second BH 
amplitude, entering the denominator, increases for ^ —\v\- Forming the moments with respect 
to the simpi (the right panel of Fig. |H)), the asymmetry is caused by the interference with the 
second BH amplitude. As explained above, we expect that this amplitude is in general larger than 
the first one and results into a rather sizeable asymmetry not only in the time- but also in the 
spacelike region. 
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Figure 8: Beam spin asymmetries of and ^lu'^'^ functions of ^, respectively, for = —0.2, 
A2 = -0.2 GeV^ y = 0.5 and different GPD models, specified in Fig.[7| 

5.5 Charge and angular asymmetries 

Let us now comment on charge and angular asymmetries, in which the Fourier coefficients cCq^^, 
cc\q, cc\2 — ss\2, CC21 — SS21, and CC32 — SS32 as well as the complementary set of the second 
interference term are attainable. As above, the Fourier coefficients of both interference terms are 
related by 
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(162) 



There are now three independent Compton form factors Cv,unp(-^)) Cv,unp{^L), and C^,unp(-^)5 
given in Eq. ()118p . which can be accessed there. Consequently, there exist two constraints among 
five nontrivial Fourier coefficients. Provided the Callan-Gross relation is assumed to be fulfilled, 
this number increases to three: 
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We add that charge and angular asymmetries can be combined with double spin-flip experiments, 
which offer information on a new combination of GPDs. As in the case of the beam-spin asymme- 
try, the number of independent Compton form factors is, however, reduced to two, as a consequence 
of our results (TOUl and (TOSl) . 
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The charge odd part is given by the interference of the first BH amphtude with the VCS as 
well as with the second BH ones. For unpolarized settings the charge asymmetry reads 

da+ - da- oc Tb^h.Tbh, + 3?e {T*j,Jycs) ■ (164) 

Taking now moments with respect to the solid angle of the final state that are even under refiection, 
e.g., by means of the weight function 

w^''^''\(j)i,9i) = {1, cos(j)iCos9i, cos{2(j)i), sin0£cos6'£, sin(20£), ...}, (165) 

the contamination of the BH interferences drops out: 

I d/2,«;^^^°(0,, Oe)^^^-^^ oc J d/2,^^^-((/),, 6^) 3?e {T*^Jycs) ■ (166) 

Corresponding to the choice of the weight function, this average will provide Fourier series in 0, 
where the zeroth, first, second and third harmonics lead to access to all leading- twist coefficients 
of the first interference term. In case when only the lepton beam of a specified single charge is 
available, on can form asymmetries with an odd weight 

^°'"(</'^,^^) (167) 
= {cos^^, cosifie, cos(2(^^) cos^^, cos{3ipi), simp^, sin(2(^^) cos sin(3(^^), ...}, 

so that the squared amplitudes exactly drop out 

J dnew°'"'{cl>t,9e)-^^cx J dQ,w'^''\(t>,,e,) {±T*^,TBn, + {T*^J^cs)} ■ (168) 

After the subtraction of the remaining BH interference is done, one measures the leading twist-two 
Fourier coefficients. Still, this procedure may allow a handle on the real part of the Compton form 
factors. If both kinds of the lepton-beam charges are available, the BH contribution drops in the 
charge even combination 

j df2ew°''\<j)e,ee)^^^^^^ J dQ,w''''\(t^,,e,)^e{r*^J^c^) . (169) 

To illustrate the feasibihty of the subtraction procedure, we consider the charge and the angular 
asymmetries 

l^cosv^, J AT 7^/4 Jo Jo ^ ^\ da-/d(2ed(p J' ^ ^ 
performed with respect to 2 cos (pi, where in both cases we choose the normalization to be 

/•37r/4 /•27r /•27r (]fj- 

/ dOt / dcp / d^t-TTT-n- 
Jn/4: Jo Jo aii(_a(p 
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Figure 9: The charge asymmetry A^^^'^ and angular asymmetry A'^°^'^'^ are displayed versus ^ 
in the left and right panel, respectively, for the same kinematics as in Fig. |H1 and different GPD 
models, specified in Fig. [7| 

These asjnumetries project the Fourier coefficient cc^q of the second BH-VCS interference term, 
which is, in the absence of J-'l, proportional to 

cc% oc 3fte I ^FiH - --^F2£ + r/ (Fi + F2) h] , (171) 

making use of Eq. ()B.28|) . One realizes that 7i is now suppressed by a factor ^/rj and, thus, with 

decreasing |^| the contribution of Ti starts to be important. However, for the sake of simplicity, 

we will set Ti for the illustration purposes to zero. In Fig. El we display these asymmetries in two 

different panels for the same kinematics and GPD models as employed above for spin asymmetries. 

Clearly, the shape of these asymmetries is in any case dictated by the GPD models. The left 

panel shows the charge asymmetry. One realizes, that as in Fig. ((7j), the D-term is responsible for 

the sign change of the asymmetries^, when C, approaches ^1]: compare solid (dash-dotted) with 

dashed (dotted) line. The size of the asymmetries in the range —\ri\ < ^ < \ri\ is driven by the 

parameterization of the reduced GPD, see also Fig. [7| We observe that in the timelike region 

(^ < 0) both asymmetries are rather similar, while in the spacelike region the angular asymmetry, 

compared to the charge asymmetry, is shifted downwards. This is caused by the interference of 

both BH amplitudes, where the first one becomes small, compared to the VCS amplitude, in the 

timelike region. Due to the additional power of ^, indicated in Eq. (jl71|) . the charge asymmetry 

with the DD-model goes to zero at ,^ ^ 0. This features will not be changed, if one includes Ti, 

since its real part is antisymmetric in ^. For the FPD-model both asymmetries go to zero which 

is caused by the fact that the normalization A/" as in the case of beam-spin asymmetries strongly 

^We remind the known fact, that for ^ ~ 1] there is a competition in sign between the regular valence and sea 
quark GPDs. The resulting sign at ^ = ±77 is a consequence of chosen parametrizations, and can not be taken as 
a clear-cut signature for the D-term contribution. 
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increase for ,^ 0. We note that the wriggles around the point ,^ ^ of the sohd and dashed 
hne in the left panel arise from the competition in the numerical increase of the numerator and 
denominator. 

We remark that as in the case of DVCS charge and angular asymmetries might be contaminated 
stronger by twist-three effects than the beam-spin asymmetries. They are mainly of kinematical 
origin, i.e., expressed by twist-two GPDs and generate Fourier coefficients cc™^, ss™""" ~ cc™""", 
and ss™""" ~ cc™""" that do not necessarily vanish at the kinematical boundaries. 



5.6 Collider experiments 

Finally, we address the question of whether GPDs can also be measured in the small-i] region 
with lepton-hadron collider experiments. Analogously to the above Eqs. (jl53p and (jl54|) one can 
estimate the ratio of the BH and VCS cross sections by 

A{^ + r])il-y)-y\^~r])smH, Fi(A2)2(l + . . .) " 

Here the ellipses stand for the ratio of terms that arise solely from the scattering by longitudinally 
and transversely polarized photons. They vanish in the DVCS limits and will be omitted since 
they do not affect our qualitative considerations. The first kinematical factor, essentially equal 
to |A^|/(Q^ + M.1i)i should be much smaller than one to ensure the applicability of perturbative 
QCD. The second term is, over a wide kinematical range, (much) smaller than one, except when 
one approaches the spacelike DVCS limit. In this case it goes like l/y^ and gives the desired 
kinematical enhancement, which allows to measure the DVCS cross section at small value of y. 
Staying away from this limit we already see that the VCS cross section is kinematically suppressed 
by at least one order of magnitude, or so, compared to the BH one. Of course, the behavior of 
rjT-C{^ , 7] , A"^) is completely unknown in the central region, since it has never been measured so 
far. Theoretically, we are facing here new non-perturbative phenomena for which the connection 
to forward parton densities is lost. So there is a lot of space for speculations, which should be 
clarified by experimental measurements. If we take the FPD-model, then 7i(^, rj, A^) grows like 
^^^^^ as C, goes to zero. Certainly, the kinematical suppression will be overwhelmed by this rise 
in the equal mass limit, ~ M^^. We have found no plausible arguments of why this should 
happen. We note, that if the sea-quark GPD does not vanish at ^ = as x = 0, the real part of the 
Compton form factor will blow up. In case the GPD is vanishing at this particular point, we still 
expect that the real part, given by the integral ()157p over the whole x range, goes with 
which follows from rather general considerations based on the DD-representation of GPDs and 
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the fact that they have to reduce to the parton densities in the kinematical forward hmit. The 
available DVCS measurements, analyzed with leading order formulas of perturbative QCD, tell us 
that the sea-quark GPD is rather sizeable at the point ^ = —r]. Consequently, there is only a small 
amount of the phase-space left, in which the GPD goes to zero at x = 0. The way this happens 
will be essentially up to the prefactor we do not know. We conclude, that the only chance to 
have an experimental handle on the VCS cross section is in the kinematical region in which 
is small, + A^£^ should be several times larger than the average ((— A^)) and at least one or 
two GeV^ and y should rather large to ensure that |?7| is small. For instance, in the case of HI 
and ZEUS experiments we find for y = 0.5, r] = -2.2 x 10-^ ((A^)) = -0.1 GeV^ and 9e = n/i 
within the DD model that the ratio ()172|) is of order one for a large interval of —\r]\/2 < ^ < \ri\. 
Similar estimates can be done for the charge and angular asymmetries, where one should project 
the Fourier coefficient Cq^, which, compared to cf^, is not plagued from an additional suppression 
factor C^/t]. 

6 Conclusions 

In the present paper we have studied the process eN — > e'N'ii. We have elaborated the structure 
of the cross section to leading power in the hard momentum. The power suppressed twist-three 
contributions, though have not been discussed presently, generate further harmonics in the squared 
amplitude. For instance, they induce the off-diagonal elements in the coefficient matrix of the 
squared VCS amplitude, e.g., ccnm with n ^ m etc., however, they will not contaminate already 
existing Fourier harmonics, e.g., cc„„ etc. This process, we have discussed, is the most favorable 
for experimental measurements of GPDs by a number of reasons: 

• it is a clean electromagnetic process which does not involve other unknown non-perturbative 
function and, thus, has no contamination from other unknown sources; 

• the virtuality of the final state photon allows to disentangle the dependence of GPDs on 
both scaling variables and thus constrain the angular momentum sum rule; 

• studies of the angular dependence of the recoiled nucleon and of the lepton pair are com- 
plimentary and lead to a rich angular structure of the cross section that can be used for 
separation of diverse combinations of GPDs; 

• variation of the relative magnitude of space- and timelike photon virtualities allows to access 
distributions of partons and anti-partons in the "exclusive" domain ,^ > \ri\. 
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• the higher one goes in skewness rj, whose maximal value is limited by the magnitude of the 
momentum transfer to be within the region of applicability of QCD factorization |A^| <^ p-q, 
the more surface in the exclusive domain (|?7| > |^|) one measures in experiment. This 
diminishes the uncertainty coming from the inaccessible inclusive sector (Ir^l < |^|). The 
exclusive domain might saturate the spin sum rule even for moderate t] since the second 
moment required for it is not extremely sensitive to the large-^ behavior of GPDs where 
the latter is known to decrease according to the quark-counting power law of conventional 
parton distributions. 

Another interesting feature of this process is that the zero value of generalized Bjorken variable 
can be exactly attained when the incoming and outgoing photons have about the same absolute 
values of virtualities ^ M^^. 

We have presently discussed the most favorable observables, namely, diverse lepton-spin and 
azimuthal asymmetries that are sensitive to the imaginary part of the Compton form factors and, 
thus, directly to GPDs. We have not discussed in full, however, phenomenological consequences 
of polarized targets, though, we have derived the complete set of formulas with explicit angular 
dependence which can be used to extract complimentary combinations of Compton form factors 
from experimental data. Longitudinal and transverse nucleon-spin asymmetries combined with 
the Fourier analysis will serve this purpose analogous to the lepton-spin asymmetries addressed 
presently. For the complete analysis along this line of exclusive electroproduction of a real photon 
see Ref. The process eA^ —>■ e'N'ii with both photons being virtual is unique due to the 
independence of skewness rj from the generalized Bjorken variable ^. Unfortunately, it suffers from 
very low cross sections, however, this drawback will be circumvented with future high-luminosity 
machines. The current analysis of available events from CLAS detector at Jefferson Laboratory 
is under way P^ . 

We would like to thank A.V. Radyushkin for useful discussions on the factorization issue and 
M. Diehl, M. Guidal, and M. Vanderhaeghen for discussions on phenomenological aspects of the 
paper. We thank the Department of Energy's Institute for Nuclear Theory at the University of 
Washington for its hospitality during the program "Generalized parton distributions and hard 
exclusive processes" and the Department of Energy for the partial support during the completion 
of this paper. The present work was supported by the US Department of Energy under contract 
DE-FG02-93ER40762. 
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A Light-cone vectors 



The vectors defining the kinematics of the process can be used to construct a pair of the light-cone 
vectors and n*, such that = n*"^ = and n ■ n* = 1, as follows 



= — / qt^ + — ^^^^^Pm ■ 



where 5^ = (M^ — A^/4)/Q^. Throughout the paper we use the following decomposition of a 
given four-vector f ^ in its light-cone components 

= v_n^ + w+n* + v±^ , (A.2) 

so that a scalar product is written as v ■ u = v^u^ + v_u^ — v_i_ ■ u±. 
The light-cone decomposition of the momenta reads 

i + v/i + 4(e5f 1- v/i + 4(e5f 

Setting Ax = we get from the last equation the minimal value of A^, namely, 



Amin = 1^ { 1 - + ^Mie/Q' - + ^Mie/Q'f + IQMievVQ'] ■ ( A.4) 

B Results for polarized target 

In the body of the paper we have given the results for the unpolarized nucleon target only. In 
the subsequent two appendices we fill the gap left and extend them by including polarization. We 
parametrize the polarization vector of the nucleon by a polar and an azimuthal angles and ^, 
respectively. In the target rest frame TRF-I it has the form (see Fig. 12)) 

S = (0,sin6'cos<P, sin6'sin<P,cos6'). (B.l) 

Introduction of the transverse polarization results into addition of an extra integration variable in 
the phase space given by Eqs. Q and 

d<P 

C/LIPS4 — X rfLIPS4 . (B.2) 
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The Fourier coefficients ab^^ = (cc^„, cs^^, sc^,^, ss^^) depend on the nucleon polarization vector 
()B.1|) and admit the decomposition 

^^nm = ^^nm,unp + ^os ab^^ Lp + sin ab^^ j.p(<P) . (B.3) 
B.l Squared VCS amplitude 

The general structure of products of Compton form factors reads 



^ VVt ^ CVgLp + C",1p cos + - %f C^v,TP sin {<P - <P) sin , (B.4) 



Iv^t , c^cs^^^ + c^^^^ eos + y \/ 1 - %f / ^ C^AT. cos - 0) sin , (B.5) 



with AA^ product having the decomposition similar to Eq. ()B.4|) . Here the polar and azimuthal 
$ angles parametrize the nucleon polarization vector 5* according to Eq. ()B.1|) . The nonvanishing 
C-coefficients in the squared VCS amplitude, defined in Eq. (jB.4IIB.5|) read for the polarized nucleon 
target 



c^v,Tp(-^,^*) = (1 - ^) ins* - sn*) , (B.7) 
cIS'tp(-^, n = (1 - v)v {nr - sn*) , (b.s) 

CvStp(-^, = -(1 - V) (vnS* + SW) + rfSS* . (B.9) 



We note that with this definitions, we have the following correspondence with those given in Ref. 
for polarized target in the DVCS limit r] = 

/7DVCS/-7r -77* ^ DVCS ^ VCS / -r ■Tr*\ , r^CS / -r* t\ 
'-'LP K-^^-^ J — ^VA,LP\'^ '-^ J ^ ^VA,LP\'^ ■,-')■, 

l-'TP+ I-'" ) J — '-'V^,TPI-'" 5 J "T '-'V.A.TP I''" ; J ) (^o.iUJ 

pDVCS/^ ^^*^ DYps ^VCS / -r- -r•*^ i /^VCS / -r•*^ 
L-rpp_ 1^7" , 7" J — L'VV,TPI'' 5 J "T '-'^y4,TP I*' I-' J- 

The remaining Fourier coefficients read for longitudinally polarized target 
cc^^l'p = 2A(2 - y)yi2 - 2y + f) {C^%A^, ^ + C^^pi^^*, T) } , 
cc^Sp = Ayy(2 - y)^{\ - y){\ - - v') {C^Zpi^^L^n + CSp(-^2'-^)} ' (B-H) 

and for transversely polarized target 



CL00,TP 



M2 V A2 V 1 - 
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cc 



vcs 

11, TP 



CC 



VCS 
22,TP 



X 



A cos {<P - 0) (2 - y)y{2 - 2y + f) {C^^%A^^ + C^a%A^\ ^)} 

+ I sin - 0) (2 - 2|/ + y^) (2 - 2y + y^) {C^^%A^. + CIS'tp(-^, ^*)} 



+ 82 sin {<P - (j)) {I - y){l - y) 



^ 'Ir-YCS f-rr -77* ^ 

L'VV.tpI-' l) 



X |a cos - 0) 1/(2 - {C^5;5rp(-^L, -^*) + C^ATpin. ^)} 
+1 sin - 0) (2 - |/)(2 - {C^5^TP(-^L, + C^5tp(-^I, •^)} 



A2 



1 



A? 



1 + r7 



/ M2 V A2 V 1 - r/ 
x4z sin - 0) (1 - y) (1 - y) {C^5'tp(-^, ^ ' CXS'tp(^, •^*)} 

where according to the general twist-two relation 



vcs _ vcs 

^^ll,LP/TP — '-'-ll.LP/TP ; 



vcs ^ VCS 
^^22, TP — ^^22, TP 



(B.12) 



(B.13) 



(B.14) 



(B.15) 



B.2 Interference term 

Similarly to the previous appendix, we present here the results for the interference term on a 
polarized target. 



\S2\ yyi \\yK\ 



COS0I /-rJ^i^' 

-Cv,unp(^) < ) - COSfc' zCv,Lp(^) < 

COS (jDf sin g 



sin 



, , , sm0 , , , . cosc, , 

iCv,tM^) { ) cos(^ - 0) + iCv,TP-(J^) < ^ sin(<? 

smipA [cos if A 



|7^2^ \\yK\ 



«Cv,unp (^) S ^ + COS O Cv,LP (^) < 

sin ix>t cos g 



+ sin 6* 



m 

\yL\ 



, , , COS0 , , , , smv , 

Cv,TP+(^) <! ^ cos(^ - 0) - Cv,TP-(-^) <^ y sin(^ - 

cosg^J IsmgjH 



sm I 



sm ( 



cos6'C^,Lp(^) 



COS( 



cosg9£ I 



sin 



m 

\vl\ 



cos ( 



C^,TP+(-^)<{ "^C0S(<? 
COS 



sm I 



Catp-(-^)<^ "^sin(^ 

Ism ^pi\ 



fB.16) 



(B.17) 



(B.18) 
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|7^l 
|7^2J 



A = AQ 



+ sin 6* 



^yl -ri)iq I \yK\ 



\yK\ 



C^,unp(-^) 



m 

\yL\ 



iC 



-4,TP+ 



,unp 

sin( 
I sin 



cos (p 
[COS If i 

cos(<P 



+ cos6'zC^,Lp(^) 



sm I 
sin g 



(B.19) 



+ iC 



ATP- 



f COS0 


> sin(<P - (j)) 






) 


[cos ip^^ 







where we used in analogy to the definitions the shorthand notation 



L 
L 



1 J{l-y){i-r^) 



(B.20) 



For the longitudinally polarized nucleon, we get the following combinations of the electromag- 
netic and Compton form factors: 



C^,LP = FiH - T] 



1 — rj 



fB.211 



N 



while for the transversal case we have four more combinations 



Cv,TP+ 
C^,TP+ 
Cv,TP- 
C^,TP- 



2'q 



AMI 



^^F, ( n - ) - ^ ((1 - ^^)F,n + r^(Fi - vF,)£ 



1 — 1] 

^1 

1-^1 

_2rf_ 
1 — 1] 



ri'F, {H + £) 



1 — 1] 
A' 



AMI 



AMI 



'l-if)F2'H-{Fi+ifF2)£ 



(B.22) 
(B.23) 
(B.24) 
(B.25) 



First, we display the explicit form of the Fourier coefficient in the interference term with the 
second BH process for an unpolarized nucleon target, which we deduced by symmetry considera- 
tions in the main text. They are 



cc, 



OljUnp 



CS, 



OljUnp 



8K(2 -y){l- y){2 - y)^-±l^e\cv ,nnA^) - C^,unp(^) 

9,\Ky{l -y){2- y)^^ 53m|-Cv,unp(-^) + C^,unp(-^) - 



Cv,unp(^L) kB.26) 



Cv,unp(^L) ,(B.27) 



cc 



lOjUnp 



-m 3fte I (2 - 2y + 2/2) (2 -2y + f) (-Cv,nnA^) - C^,unp(^) 



1] 

Kl - Z/)(l - y)^— ^Cv,unp(^L) , 



i]i 



cc 



sc 



12,unp 
2 

10,unp 



-1QK{1 -y){l- yP '-^e {Cv,u„p(^) - C^,unp(-^)} 



8XKy{2 - y){2 - 2y + y')^m |Cv,unp(-^) - ^C^,unp(-^) 



(B.28) 
(B.29) 
(B.30) 
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CC^l,,np = 8K{2 -y){l- m - y)^3fJe|Cv,unp(^) + C^,unp(^) - ^Cv,unp(-^L)|(B.31) 
SC^l,,np = -^\Ky{l - y)(2 - 53m|cv,unp(-^) + C^,unp(-^) + ^Cv,„np(^L)| ,(B.32) 



(B.33) 



Next we consider the Fourier coefficients for a polarized target. From the results ()B.16|) - ()B.19|) 
we immediately read off several relations which allows us to obtain the expression for the Fourier 
coefficients making use of simple substitution rules in the unpolarized case (|119p - (|126|) and ()B.26|) - 
(ESI: 



r tINT r tINT 

|cCoi,CCio,CC2i|lp = |cSoi,SCio,SC2l|, 



I unp 



SmCunp^SReCLP ' 



(B.34) 



r iINT r T 

jcSoi, SCio, CS12, SC21, SC32ILP = jcCoi, CCio, CC12, CC21, CC32I 



INT 

unp 



5RcCunp^3mCLP 

In the case of a transversely polarized target we have as before an additional decomposition in 
cos(^ — 0) and sin(^ — 0) and should replace K {K) by L (L) 



{CCOI, CCio, CC21}™ + = C0S{<P - (p) {CSOI, SCio, SC2l}™'^ 



I unp 



jcSoi, SCio, SC2i}tp_ = Sm(0 - 0) jcSoi, SCio, SC21I 



{cSol, SCio, CS12, SC21, SC32}?^>+ = C0S(^P - 0) {cCol, CCio, CC12, CC21, CC32} 



unp 



GmCunp^9?eCrpp^ , 



(B.35) 



{cCoi,CCio,CCi2,CC23,CC32} 



INT 
TP- 



sin(^P - 0) {ccoi, ccio, cci2, CC23, CC32} 



INT 
unp 

INT 

unp 



McCunp^3mCrpp_|_ , 



INT ^^INT „^INT 
21 5 



Applying the relation (jll3p . the remaining nonvanishing Fourier coefficients, 
cs™^ for LP and TP+ as well as ss™"'", ss™"'", ss™"'" and cs™""" for TP— are easily established. 
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